Deviations from Keplerian Orbits for Solar Sails 



Roman Ya. Kezerashvili and Justin F. Vazquez-Poritz 

Physics Department, New York City College of Technology, The City University of New York 

300 Jay Street, Brooklyn NY 11201, USA 

The Graduate School and University Center, The City University of New York 
365 Fifth Avenue, New York NY 10016, USA 

rkezerashvili@citytech . cuny . edu jvazquez-poritz@citytech . cuny . edu 

o 

O i ABSTRACT 

It is shown that the curvature of spacetime, a possible net electric charge on the sun, a small positive cos- 
mological constant and the oblateness of the sun, in conjunction with solar radiation pressure (SPR), affect 
the bound orbital motion of solar sails and lead to deviations from Kepler's third law for heliocentric and 
non-Keplerian orbits. With regards to the Lense-Thirring effect, the SRP increases the amount of precession 
per orbit for polar orbits. Non-Keplerian polar orbits exhibit an analog of the Lense-Thirring effect in which 
the orbital plane precesses around the sun. 
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INTRODUCTION 

The widely-discussed concept of the Solar Probe mission is an integral part of the programs of many 
national space agencies. No previous spacecraft has traversed the region of less than 0.3 AU from the sun. 
However, during the last two decades, a number of solar probes have been designed to perform in-situ near- 
sun measurements in a rather unusual and difficult environment. With the success of the Ulysses mission and 
the unprecedented opportunities offered by the SONO mission launched in 1995, a golden era has dawned in 



solar physics. It is therefore not surprising that a solar probe penetrating to a few solar radii has once again 
come into focus [1]. The primary goals of the solar probe mission are to explore the solar corona, unveil the 
formation and origin of the solar wind and understand the mechanisms for storing and transporting energetic 
charged particles. The solar probe will be a voyage of discovery to the center of our solar system. 

The primary focus of the solar probe will be to address outstanding questions regarding the structure, 
dynamics and physical properties of the outer solar atmosphere. This region remains the last unknown link 
between the solar output and the terrestrial environment. A close-solar encounter concurrently provides the 
opportunity to conduct other scientific investigations, such as the effects of general relativity on the trajectory 
of a solar sail spacecraft. In contrast with the conventional solar probe mission, a solar sail propelled (SSP) 
satellite allows for non-Keplerian orbits, for which the orbital plane does not pass through the center of mass 
of the sun and the SSP satellite is levitated above the sun. General relativistic effects on a levitated near-sun 
solar sail is of particular interest. 

One of the most basic laws that describes motion in the solar system is Kepler's third law, which can be 
derived from Newton's law of gravitation and provides a relationship between the period T and the orbital 
radius r for objects in circular orbit around the sun: 

n 47I 2 -3 

T 2 = r 3 , (1) 

GM 

where G is the gravitational constant and M is the mass of the sun. For elliptical orbits, we replace r by 
the length a of the semi-major axis of the ellipse. We will discuss deviations from Keplerian orbits due 



to phenomena within conventional physics which can be observed from the motion of an SSP satellite. In 
the general relativistic framework, objects tend to follow geodesies on curved spacetime. At the same time, 
the sun emits electromagnetic radiation which produces an external force on objects via the solar radiation 
pressure (SRP). Therefore, we can say that objects move in the photo-gravitational field of the sun. It is of 
interest to analyze how the trajectories of objects deviate from geodesies due to the SRP. We will assume 
that the backreaction of the solar radiation on spacetime is negligible. The purpose of this paper is to point 
out various sources of deviations from Keplerian orbits and to discuss how the resulting change in period is 
enhanced by the SRP to the degree in which it may be observed for some cases. 

1. GENERAL RELATIVISTIC EFFECTS FOR HELIOCENTRIC ORBITS 

We first consider the effect of the SRP on Keplerian orbits in the Newtonian approximation of gravity. 
For simplicity, we start with the scenario in which the solar sail is in circular orbit within the plane of the sun 
and the surface of the sail is directly facing the sun. Then both forces act along the same line and fall off with 
the heliocentric distance as 1/r 2 , where the SRP force is repulsive and the gravitational force is attractive. 
In this case, the mass of the sun in Kepler's third law (Q} is effectively renormalized as M ^ M = M — k/G, 
where K = i3j^=, Ls = 3.842 x 10 26 W is the solar luminosity, and o is the mass per area of the sail- a 
key design parameter for solar sails [2]. The parameter r\ lies in the range 0.5 < T| < 1, where r| = 0.5 
corresponds to the total absoiption of photons by the sail and r) = 1 corresponds to total reflection. We will 
consider an SSP satellite with the following specifications: 

r = 0.05 AUw 7.48 x 10 9 m, 

r| = 0.85, a = 0.00131 kg/m 2 . (2) 

We consider the simultaneous effects of the SRP and various general relativistic effects associated with 
the slowly rotating sun which has a mass M, angular momentum J, and possibly a small amount of net 
charge Q. The exterior spacetime in the vicinity of the sun is described by the large-distance limit of the 
Kerr-Newman metric: 



ds 2 = -fc 2 dt 2 -^ sin 2 Qdtd<\> + — + r 2 (c?e 2 + sin 2 e^ 2 ), f=l ^ + — r^- (3) 

c L r j c l r c'r 1 

where the Coulomb constant k e = 8.988 x 10 9 N m 2 /C 2 . We first consider the case of Q = 0, for which Q 

reduces to the Kerr metric. The constants of motion E = —p t /m, L z = p§/m and p 6 = jz\Jp — -J^r§ are 
associated with the energy per mass, the component of angular momentum per mass that is normal to the 
equatorial plane, and the total angular momentum per mass, respectively. P = corresponds to motion in 
the equatorial plane, for which 8 = %/2 and p® = 0. No other orbits can lie within a fixed plane. In the 
presence of the SRP, the orbital equation is [3] 

(dry_E 2 j^GM f(Lj + P) AGJEL Z (4) 
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From (O, we find that circular orbits (j£ = ^ = \ in the equatorial plane have an orbital period T of 
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We have presented T 2 as a product of three factors. The first factor is the same as in the Newtonian approx- 
imation for gravity in the presence of the SRP, and embodies the fact that the SRP effectively renormalizes 



the solar mass. The second factor shows the deviation from Keplerian orbits due to the simultaneous effects 
of the SRP and the curvature of spacetime outside of a static central body. The third factor gives the devi- 
ation from Keplerian orbits from the combined effects of the SRP and frame dragging due to the rotation 
of the sun. Note that the second and third factors in © involve both the solar mass M and the parameter K 
separately, rather than simply in the combination of the renormalized solar mass M. For a static spacetime 
(J = 0) and for no SRP effects (k = 0), © reduces respectively to 
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The first expression in ((6]) shows that it is only the simultaneous effects of the static curvature of spacetime 
and the SRP that lead to a radial dependent (and thus energy dependent) deviation from Kepler's third law. 
For the specification given in (0, we find that this yields an increase in the period of about 0.6 s. The second 
expression in ((6]) shows that, in the absence of the SRP, there is still a deviation from Keplerian orbits due to 
frame dragging. However, the SRP enhances the effect of frame dragging, since it is the renormalized mass 
M that appears in the last factor. The speed of the outer layer of the sun at its equator is v « 2000 m/s at the 
equatorial radius R rj 7 x 10 8 m, and we will assume that the core of the sun rotates with the same angular 
speed. Therefore, J = ^MvR 10 42 kg m 2 /s. Without the effects of the SRP, frame dragging leads to an 
increase (decrease) in the period of only 4 x 10~ 5 s for a prograde (retrograde) orbit. With the SRP and the 
specifications in (O, frame dragging leads to a change in period of about 0.01 s. 

The orbital plane of a polar orbit which passes through the poles of the sun will precess in the (j) direction, 
which is a well-known effect of frame dragging called the Lense-Thirring effect. Following [4], we can find 
the approximate angle of precession for a polar orbit during one orbital period up to linear order in J: 
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which is increased by the SRP due to the corresponding increase in period. For both a standard and an SSP 
satellite in a polar orbit at r = 0.05 AU, we find the rate of precession to be about 0.03 arcseconds per year. 

We will now consider the effect that a small amount of net charge Q on the sun would have on an SSP 
satellite with charge q. It has been suggested that the sun has a net charge of up to Q rs 77 C [5]. Taking 
into account the electric force, we find that 
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where we have kept only terms up to order Q 2 . Note the Q 2 term in the period, due to the backreaction from 
the sun's charge on the geometry, is present even for a neutral SSP satellite. However, this term increases the 
period by an amount of only 10~ 35 s (10~ 38 s without the SRP), which reflects the fact that the backreaction 
of the charge on the geometry is negligible, as to be expected. 

Primarily due to the photoelectric effect, but also the Compton effect and electron-positron pair pro- 
duction, a solar sail made from Beryllium, for example, will equilibrate to a charge per area of 0.065 
C/m 2 [6]. For an SSP satellite with the specifications given in [3,4] and a mass of 1000 kg, this gives 
a charge q = 5 x 10 4 C and an increase in period of about 230 s (0.05 s without the SRP). Thus, due to 
the SRP, even a small charge Q could certainly yield a measurable increase in the period, making this a 
potentially powerful test for net charge on the sun. 

Since the SRP enhances a variety of small effects, including those associated with spacetime curvature, 
one could ask whether an SSP satellite could be used to test for the presence of a cosmological constant. 
In fact, supernovae observations suggest that our universe might have a very small positive cosmological 
constant [7, 8]. In the presence of a cosmological constant A and neglecting the effects of the sun's rotation 



as well as possible net charge, the spacetime in the vicinity of the sun is described by a metric of the form 
(O, where the function / is now given by / = 1 — 2 S^- _ xhis yields 

o 47i 2 r 3 r c 2 A ,1 

T 2 w =- 1 + ^r 3 . (9) 
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For A 10~ 52 m~ 2 along with the specifications in (0, this leads to an insignificant increase in the orbital 
period of 10" 17 s (10~ 21 s without the SRP). While one might be able specify values of a and r such that 
the collection factor is actually significant, the period would then be too large to make this a feasible test for 
the presence of a cosmological constant. 

2. GENERAL RELATIVISTIC EFFECTS FOR ORBITS OUT OF THE PLANE OF THE SUN 

We will briefly review non-Keplerian orbits within the framework of Newtonian gravity. The plane of 
a non-Keplerian orbit does not pass through the center of mass of the sun, and the SSP satellite is levitated 
above the sun, as shown in Figure [TJ We orient our spherical coordinate system so that the SSP satellite 
orbits at constant heliocentric distance r and polar angle 6. A non-Keplerian orbit can be maintained with 



a SRP 




Fi gure 1 ! Orienting a solar sail at an appropriate pitch angle \|/ renders it possible to have a circular orbit outside of the plane of 
the sun, as shown in the left plot. The Lense-Thirring effect for non-Keplerian polar orbits is shown in the right plot. 

a suitable pitch angle \|/ relative to the radially outgoing solar radiation hitting the surface of the sail. This 
enables one to control the direction of thrust due to the reflected portion of sunlight. By changing the pitch 
angle, one can also transfer the SSP satellite between different non-Keplerian orbits whose orbital planes 
have different orientations, as well as between non-Keplerian orbits and orbits within the plane of the sun. 

Recall that r\ parameterizes the portion of light that is reflected by the sail. The reflected light pushes the 
solar sail at an angle \|/ relative to the radial direction, while the absorbed light pushes the solar sail radially 
outwards [9, 10]. In the Newtonian approximation for circular non-Keplerian orbits [4], 
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Now we consider the effects of curved spacetime on non-Keplerian orbits. Keeping only the leading 
correction due to spacetime curvature for orbits parallel to the equatorial plane of the sun, the period can be 
found to be 



-^sin 2 6 r 3 



GM* 



1 + 



G(M-M*} 



c 2 r 



1 + 



2-s/G J sinQ 
c 2 \[W f r i l 2 



(11) 



In dTTT > T 2 is a product of three factors. The first factor is the same as for non-Keplerian orbits in the 
Newtonian approximation for gravity. The second factor is due to the simultaneous effects of the SRP and 
the static curvature of spacetime. The third factor embodies the combined effects of the SRP and frame 
dragging due to the rotation of the sun. As was the case for orbits within the plane of the sun, the SRP 
increases the effect that frame dragging has on the orbital period, since it is the renormalized mass M* that 
appears in the last factor. 



There is an analog of the Lense-Thirring effect for non-Keplerian orbits. Namely, frame dragging causes 
the plane of non-Keplerian orbits parallel to polar orbits to precess around the sun, as shown in the right 
diagram of Figure Q] Although we do not present the resulting equations in their entirety since they are 
rather cumbersome, the angle of precession can be approximated by 

47tG/sin8 ^ 
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3. THE OBLATENESS OF THE SUN 

We now consider the deviation from Kepler's third law on non-Keplerian circular orbits due to the 
oblateness of the sun. The oblateness of the sun has recently been measured with unprecedented precision 
to be as much as J 2 ~ 9 x 10~ 6 during active phases of the solar cycle [11]. Working in Newtonian gravity, 
the external gravitational potential of an oblate spheroid is given by [12] 

GM GM " /R\n , 

v = + — £/„ - p n cose , (13) 

r r n f2 KrJ 

where J n are the multipole mass moments, R is the equatorial radius of the sun and P n are the Legendre 
polynomials. The effective mass in the first term is renormalized by the solar radiation pressure, whereas 
the multipole mass moments are not affected. A consideration of the first sub-leading term in ] 2 due to the 
oblateness of the sun leads to the relation 

, 47l 2 r 3 sin 2 6 / 3MJ 2 R 2 , , A 

T 2 ^ = 1+ _ (l-3cos 2 9 . (14) 
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Interestingly enough, the oblateness of the sun increases the period for non-Keplerian orbits within the 
angular range of approximately 55° < 9 < 90°, whereas the period is decreased for 6 < 55°. The case for 
which the circular orbit of the SSP satelite is confined to the equatorial plane can be obtained from (fT4l by 
putting 6 = 90°. As an example, without the SRP, the oblateness of the sun increases the period of orbits 
within the plane of the sun by about 0.02 s. With the SRP and the specifications given in (O, the sun's 
oblateness increases the period by about 105 s. 

For non-circular orbits, the perihelion shift per period is given by 

3kMJ 2 R 2 

M~fi (15) 

The characteristics of the precession due to the oblateness of the sun are rather different than general rel- 
ativistic precession. For instance, precession due to oblateness occurs in the reverse direction relative to 
the orbit. For a conventional satellite at r = 0.05 AU, the precession rate is about 14 arcseconds per year. 
For our SPP satellite, the precession rate increases by a factor of y/M/M to about 235 arcseconds per year. 
Although the SRP increases the rate of precession due to the oblateness of the sun, it can be shown that 
the precession rate due to spacetime curvature is decreased by the SRP. Perhaps the SRP could be used to 
disentangle these two types of precession. 



CONCLUSIONS 

We have considered the effects of various phenomena along with the SRP on bound orbits of a solar 
sail propelled satellite. The consideration of the SRP on its own leads to a renormalization of the effective 
solar mass. However, when the SRP is coupled to other effects such as spacetime curvature, the orbital 
characteristics depend on the solar mass M and the solar sail parameter K separately, and not simply on 
the combination of a renormalized mass. The SRP affects the period of the SSP satellite in two ways: by 
effectively decreasing the solar mass, thereby increasing the period, and by enhancing the effects of other 



phenomena by three orders of magnitude or more, rendering some of them detectable. When the SRP 
is not taken into account, the static curvature of spacetime does not alter the period of circular orbits as 
given by Newtonian gravity However, when the effects of spacetime curvature and the SRP are considered 
simultaneously for an SSP satellite, there is a deviation from Kepler's third law. 

There are two sources of perihelion shift for non-circular orbits. Spacetime curvature causes the perihe- 
lion to be shifted in the direction of the orbit, whereas the oblateness of the sun causes the perihelion to be 
shifted in the reverse direction. Without the SRP, the perihelion shift due to spacetime curvature is generally 
larger than that due to the oblateness of the sun. However, this is not the case when the effects of the SRP 
are taken into account on an SSP satellite. 

For orbits within the plane of the sun as well as non-Keplerian orbits, we have considered the Lense- 
Thirring effect, which is the precession of polar orbits due to the rotation of the sun. The rate of precession 
is not altered by the SRP up to leading order in the angular momentum of the sun. In particular, any satellite 
in a polar orbit at a distance of 0.05 AU from the sun will experience a rate of precession of about 0.03 
arcseconds per year. 

Observations of the trajectory of a near-solar mission would provide an interesting confirmation of these 
various phenomena. 
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